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Abstract 

The spherically symmetric Einstein- Vlasov system in Schwarzschild 
coordinates (i.e. polar slicing and areal radial coordinate) is consid- 
ered. An improved continuation criterion for global existence of classi- 
cal solutions is given. Two other types of criteria which prevent finite 
I/"") ■ time blow-up are also given. 

o 
in 

O ■ 1 Introduction 

o 

A central issue in general relativity is the problem of spherically symmetric 
£h . gravitational collapse and to determine whether or not the weak- and strong 

cosmic censorship conjectures (CCC) hold true in this case. These conjec- 
tures were raised by Penrose in the late '60s. By now they have precise 
^ ■ mathematical formulations |S] and are considered to be among the most im- 

portant problems in classical general relativity. In brief the meaning of the 
weak CCC is as follows: The Penrose singularity theorem says that there are 
initial data that will lead to singular spacetimes, i.e. geodesically incomplete 
spacetimes. The weak cosmic censorship states that generically, singularities 
will always be covered by an event horizon and distant observers will never 
be able to "see" the singularities. Singularities that can be "seen" are called 
naked. To answer the CCC a detailed picture of the global properties of 
the solutions to the Cauchy problem for the Einstein matter equations are 
necessary. 

In 1999 Demetrios Christodoulou finished his long term project on the 
spherically symmetric Einstein-Scalar Field system (ESF) by proving both 
the weak- and strong cosmic censorship conjecture jB] by a careful analy- 
sis of the Cauchy problem. This important work is the only case where 
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such a study has been completed in an affirmative sense. On the contrary, 
Christodoulou has shown that when the matter model is dust (i.e. a pres- 
sureless perfect fluid) then there are naked singularities generically |EL It is 
however a rather common belief that a dust model is not an appropriate mat- 
ter model for studying cosmic censorship since the pressure is assumed to be 
zero. A study of a more general fluid model for gravitational collapse suffers 
at present from a lack of satisfying mathematical results of such equations 
in the absence of gravity and in the case of Newtonian gravity. A successful 
treatment in these cases seems necessary before it will be mathematically 
meaningful to couple such a model to the Einstein equations. A matter 
model that is known to have satisfying mathematical features is the Vlasov 
model which gives a kinetic description of the matter. In particular it was 
shown independently by Pfaffelmoser [2j and Lions and Perthame |17j that 
in the case of Newtonian gravity, where a kinetic description gives rise to 
the Vlasov-Poisson system, classical solutions exist globally in time. This 
system can now be said to be well-understood mathematically and it is 
therefore natural to choose the Vlasov model for investigating gravitational 
collapse in general relativity. In particular this matter model is phenomeno- 
logical in contrast to a scalar field which is said to be field theoretical. It 
would certainly be desirable to establish similar results for a phenomenolog- 
ical model as Christodoulou has done for a scalar field. One way to attack 
the cosmic censorship conjecture is to prove global existence in a singularity 
avoiding coordinate system. It is then often conjectured that cosmic censor- 
ship would follow straightforwardly. The first attempt by Christodoulou on 
the ESF system was indeed of this nature. However, in that case it turned 
out that such an approach was doomed to failure |HJ by finding initial data 
that lead to naked singularities In the case of Vlasov matter it is still 
an open question whether or not naked singularities form for any initial 
data. (Note that cosmic censorship holds true if the set of initial data which 
generates naked singularities is a null set). Thus the approach of proving 
global existence in a certain coordinate system is yet a possile path to weak 
CCC. The aim of this paper is to present some new results on the problem 
of global existence of the spherically symmetric Einstein- Vlasov system in 
Schwarzschild coordinates. In particular we give an improved continuation 
criterion. We also consider two other types of criteria, one in which we study 
the case when all matter is assumed to fall inwards. A continuation criterion 
for this system was first given by Rein and Rendall in [2Sj. To discuss the 
content of this criterion we need to introduce the main concept in kinetic 
theory, the distribution function. Consider a collection of particles, where 
the particles are stars, galaxies or even clusters of galaxies. A characteristic 
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feature of kinetic theory is that the model is statistical and that the particle 
system is described by a distribution function / = f(t,x,p), which repre- 
sents the density of particles with given space-time position (t, x) £ R x M 3 
and momentum p£R 3 . Let us define the function Q(t) by 

Q(t) := sup{|p| : 3(s,x) G [0,t] x M 3 such that f(s,x,p) / 0}. (1) 

Hence, there are no particles at time t having larger momentum than Q(t) 
in the system. If we now consider a solution of the spherically symmetric 
Einstein- Vlasov system on a time interval [0, T[, then it is proved in [23] 
that if Q(t) is bounded on [0, T[ then the solution can be extended beyond 
T. This is called a continuation criterion. This is analogous to the situation 
for the Vlasov-Maxwell system where the same condition ensures that the 
solution can be extended. This was proved in a classical work by Glassey 
and Strauss ^H] but recent results by Pallard |20j improve this criterion. 

The outline of the paper is as follows. The Einstein- Vlasov system in 
Schwarzschild coordinates is introduced in section 2. A brief review of pre- 
vious results is also given here. In section 3 we present some general a priori 
bounds not previously in the literature. In section 4 we make a compari- 
son between the Einstein- Vlasov- and the Vlasov-Poisson system to point 
out the main differencies. An improved continuation criterion is proved in 
section 5 using a mixed L 3+ and weighted L°° approach. In section 6 an 
approach based on hypotheses along characteristics is given. A discussion 
about what is expected to happen if these bounds are not satisfied is also 
included. In section 7 the case when all matter is falling inwards is studied. 
It is seen that the well-known result that a Schwarzschild spacetime with 
mass M has closed null rays when r = 3M, also plays an essential role for 
the Einstein- Vlasov system. 

2 The Einstein- Vlasov system 

For a derivation of the system given below we refer to |22j and |25j . In 
Schwarzschild coordinates the spherically symmetric metric takes the form 

ds 2 = - e Mt,r) dt 2 + e 2X(t,r) dr 2 + ^2 + ^2 Q^y (2) 

The Einstein equations are 

e- 2X (2r\ r - 1) + 1 = 8irr 2 p, (3) 
e~ 2x (2rnr + !)-! = 8irr 2 p, (4) 
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X t = -4irre x+ »j, (5) 

e" 2A (/irr + (jir - K)(Vr + -)) - e"^(A tt + A t (At - //*)) = 4vrg. (6) 

r 

The Vlasov equation for the density distribution function / = f(t,r,w,F) 
reads 

d t f + e^^drf - (X t w + e^V^ - e^ x -S^)d w f = 0, (7) 

where 



£ = E(r, w, F) = a/1 + w 2 + F/r 2 . (8) 
The matter quantities are defined by 



oo />oo 



p(i,r) = 4/ / Ef(t,r,w,F) dwdF, (9) 

r J-ooJO 



P(^ r ) = "a/ / -jjrf(t,r,w,F) dwdF, (10) 

i/ — oo «/ 

/'OO /'OO 

j(t,r) = - / / wf{t,r,w,F), dwdF (11) 

r J-ooJO 

/*oo /-OO ^ 

g(i,r) = -4 / / -f(t,r,w,F)dwdF (12) 

i/ — oo «/ 

For technical reasons we also introduce the quantity 

f'OO f'OO 

k(t,r) = — / H/(t,r,«;,F) dwdF, (13) 

r ./-oo J 

so that k is nonnegative and dominates p and \j\ but is dominated by p. 
The variables ui G (—00,00) and F G [0, 00) are the radial component of 
the momentum and the square of the angular momentum respectively. The 
following boundary conditions are imposed 

lim \(t,r) = lim /j,(t,r) = A(i,0) = 0, Vt > 0. (14) 

r— »oo r— >oo 

The last condition follows by requiring a regular center. We point out that 
the Einstein equations are not independent and that e.g. the equations (JSJ) 
and © follow by ©, ® and 0. 

For an extensive review of previous results of this system see pQ. Here 
we give a short summary of the results concerning the Cauchy problem. 
The issue of global existence for this system was first investigated by Rein 
and Rendall in where they obtain global existence for small (compactly 
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supported) initial data and in addition they show that the resulting space- 
time is geodesically complete. They also prove the continuation criterion 
mentioned in the introduction. In j^Tj Rendall studies the same system 
using a different set of coordinates (maximal slicing and isotropic radial co- 
ordinate). Global existence for small initial data is obtained as well as an 
equivalent continuation criterion. In [25] Rendall shows that there exists 
initial data leading to trapped surfaces and therefore to singular spacetimes 
by the Penrose singularity theorem. The natural question is then if the 
weak CCC holds. To investigate this it is necessary to consider large initial 
data. This is done in and |2Zj • I n the former case Schwarzschild coordi- 
nates are used and it is proved that if matter stays uniformly away from the 
center of symmetry then global existence follows. In the latter case where 
the maximal and isotropic gauge is used a similar result is shown assuming 
an additional condition on one of the metric coefficients. The advantage 
being that the result is obtained in a more direct way and the analysis of 
the Vlasov equation is not necessary. Recently Dafermos and Rendall jllj 
have shown a similar result for the Einstein- Vlasov system in double-null 
coordinates which has a very interesting and important consequence. In the 
proof of the CCC for the ESF system Christodoulou analyzes the formation 
of trapped surfaces and their presence play an important role in his proof. 
Dafermos ^U] then strengthens the relation between trapped surfaces and 
weak CCC; it is shown that a single trapped surface or marginally trapped 
surface in the maximal development implies weak CCC. The theorem in jlUj 
requires certain hypotheses on the solutions of the Einstein matter system. 
The main purpose of jllj is to show that these hypotheses are satisfied for 
the spherically symmetric Einstein- Vlasov system, which is indeed shown to 
be the case. This result might thus be useful in the search for a proof of the 
weak CCC. Two numerical investigations have been carried out where the 
primarily goal was to study critical collapse and determine if Vlasov matter 
is of type I or type II, i.e. if arbitrary small black holes can form (type II) 
or if there is a mass gap (type I). In the first study [201 Schwarzschild coor- 
dinates were used and in the second ^S] maximal slicing and an areal radial 
coordinate was used. Both of these studies indicate that Vlasov matter is of 
type I, and moreover, no signs of singularity formation could be seen in ei- 
ther of these two cases which support the idea that global existence holds in 
these two coordinate systems. The investigation in this paper concerns the 
Einstein- Vlasov system in Schwarzschild coordinates but we plan to investi- 
gate other choices of gauge conditions both analytically and numerically in 
the future. 

Let us now specify the sets of initial data that will be considered. Regular 
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initial data is denned as the set of nonnegative compactly supported smooth 
functions /o = fo(r,w,F), such that 



4tt 2 / / / Ef (r,w,F)dwdFdr < -. 
JO J-oo Jo 2 

The last condition says that no trapped surfaces are present initially. Given 
R- > 0, we define the class of regular initial data with radial cut-off, Ir_ , as 
the subset of regular initial data such that /o = when r < f?_. Likewise, 
given F- > 0, the subset Ip_ such that / = when F < F_, is called regular 
initial data with angular momentum cut-off, and consequently, is the 

subset of regular initial data with radial and angular momentum cut-off. 
Thus, for initial data in Ir_ f~ there is no matter in a ball around the centre 
and there is a lower bound on the angular momentum F. By conservation of 
angular momentum the latter property is then true for all times. (Angular 
momentum cut-off has shown to be essential also in previous works, see |14j 
and jHJ.) We also make the definitions 

F + := sup{F : 3(r, w) G M 2 such that / (r, w, F) / 0}, (15) 

and 

R + := mf{R : / (r, •) = for all r > R}. (16) 

Let us write down a couple of known facts about the system ([5 ]) -(|14 )l . A 
solution to the Vlasov equation can be written 

f(t, r, w, F) = / (i?(0, t, r, w, F), W(0, t, r, w, F), F), (17) 

where R and W are solutions of the characteristic system 

^ _ J»-\)(s,R) W 

ds E{R,W,FY [ ' 



dW 

ds 



-X t (s, R)W - e^- x)(s ' R) fi r (s, R)E(R, W, F) 

(»-\)(s,R) F (iq\ 

R S E(R,W,FY { ' 



such that the characteristic (R(s, t, r, w, F), W(s, t, r, w, F), F) goes through 
the point (r, w, F) when s = t. This representation shows that / is non- 
negative for all t > 0, H/ 1 |oo = ll/olloo) an d that for regular initial data 
with angular momentum cut-off, f(t, r,w,F) = if F < F-. There are two 
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known conserved quantities for the Einstein- Vlasov system, conservation of 
the number of particles 

Air 2 e A(i ' r) ( / / f(t,r,w,F)dFdw) dr, 

Jo \J-ooJo J 

and conservation of the ADM mass 



DC 



M := 4vr / r 2 p(t,r)dr. (20) 
Jo 

These conservation laws follow from general argments but can easily be 
obtained by taking the time derivative of the integral expressions and make 
use of the Einstein equations and the Vlasov equation. The mass function 
m is defined by 

pr 

m(t,r) := 4ir / rj 1 p(t,r))drj, (21) 
Jo 

and by integrating (JIJ) we find 

g _ 2A(t>r) = 1 _ Mt,r)_ 
r 

Thus, as long as the solution exists we have 

^1 < 1/2. (23) 
r 

A fact that will be used frequently is that 

p + A < 0. 

This is easily seen by adding the equations © and (j3J, which gives 

A r + p r > 0, 

and then using the boundary conditions on A and p. From (|22(1 we see that 
A > 0, and it follows that p < 0. Finally, we note that in a local 
existence theorem is proved and it will be used below that a classical (or 
regular) solution exists on some time interval [0,T[. 



3 General a priori bounds 

In this section we show a few a priori bounds which are easy to derive but 
cannot be found in the literature. Together with the conserved quantities 
mentioned in section 2, and the inequality (|23|). these results constitute the 
only known a priori bounds for the spherically symmetric Einstein- Vlasov 
system. 
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Lemma 1 Let (/,//, A) be a regular solution to the Einstein- Vlasov system. 
Then 

POO 

/ 4tt(p + p)re 2 V +A (fr < 1, (24) 
Jo 

f'°° m 

/ (4 + 4vrp)re 2 Vffr < 1. (25) 
Jo r 

Proof. Using the boundary conditions (|14|) and that ^ + A < we get 

1 > l- e ^ +A (t,0) = / ^e» +x dr 

Jo dr 

poo 

= / (/i r + A r )e^ +A dr. 
JO 

The right hand side equals (|24|) by equations © and (@J) which completes 
the first part of the lemma. The second part follows by studying instead 
of e^ +A . 

□ 

Next we show that not only p(t, •) G L , which follows from the conservation 
of the ADM mass, but that also e 2X p(t, •) G L 1 . 

Lemma 2 Let (/, p, A) be a regular solution to the Einstein- Vlasov system. 
Then 

POO POO £ 

/ r 2 e 2X p(t,r)dr < / r 2 e 2X p(0, r)dr -\ . (26) 

Jo Jo 87r 

Proof. Using the Vlasov equation we obtain 

8t (r 2 e 2X p(t,r)^J = -0 r (rV +A j) - re^ +A 2je 2A 



rn 
r 



< -8 r (rV +A j) + ^re^ +x (p + p)e 2A . 

Here we used that m/r < 1/2 together with the elementary inequality 2\j\ < 
p + p, which follows from the expressions In view of (|18|) we see 

that lim r ^ oc r 2 j(t, r) = 0, since the initial data has compact support. Since 
the solution is regular and hence bounded the boundary term at r = 
also vanishes (as a matter of fact spherical symmetry even implies that 
j(t,0) = 0). Thus, by lemma 1 we get 

poo 2 
d t / r 2 e 2X p(t,r)dr < — , 
Jo 8?r 

which completes the proof of the lemma. 
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□ 



A natural question is if the inequalies l|24|) and (|25|) can bound any matter 
quantity in an LP— space with p > 1, in view of the fact that merely r and 
not r 2 is present in the integrand. For the Vlasov-Poisson- and the Vlasov- 
Maxwell system such LP estimates (with p > 1) play an essential role in the 
global existence proofs (although the problem in 1 + 3 dimensions is open for 
the Vlasov-Maxwell system affirmative results exist in lower dimensions |13j ) 
and such results are thus desirable. We have the following result (where 
L p {R ?J ,ijj) means the LP space with weight u>): 

Lemma 3 Let (/, /x, A) be a regular solution to the Einstein- Vlasov system 
for initial data with angular momentum cut-off. Then 



Proof. Here we use the momentum variable v instead of (w,F). These are 
related as follows 



g(V) e L 5 / 4 (i? 3 , e fl+x e 2X ). 



x ■ v 



F = x 2 v 2 — (x ■ v) 2 = \x x v 



w = 



r 



where r = From it follows that q now takes the form 




By the assumption on the initial data we have that r\f\ 
support of / and we get 





(27) 



(28) 



(29) 



(30) 



(31) 
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Here we took 

K - ( f fM + Ej"" 



IRS 



r 



Now since F < F + we get that 

g^<c[ fVTT ^ dv < C P ~, (32) 
and from 1)24(1 we obtain 

poo roo 

/ r 2 i^e^ x e 2X dr <C pre 2X e^ +x dr < C, 
Jo Jo 

and the lemma follows. 

□ 

Remark. The first factor e^ +A in the weight e /i+A e 2A is the square root 
of the determinant of the metric and can thus be viewed as a part of the 
integration measure. The second factor e 2A is always greater than one which 
then sharpens the bound for q, in analogy with lemma 2. 

4 A comparison between the EV- and the VP sys- 
tem 

It is instructive to compare the structure of the spherically symmetric Vlasov- 
Poisson system, where the gravitational field is Newtonian, with the spher- 
ically symmetric Einstein- Vlasov system in Schwarzshild coordinates (in 
other coordinates such a comparison can be slightly different). The for- 
mer system is the limit of the latter as the speed of light goes to infinity, 
as was proved in [51], but the mathematical understanding of these two 
systems are nevertheless very different. In this section we try to point out 
the main differencies between these systems from a mathematical point of 
view. A Newtonian gravitational system is modelled in kinetic theory by 
the Vlasov-Poisson system which reads 

d t f + p-V x f-E(t,x)-V p f = 0, (33) 
£ = V0, A<P = p:= J Bjj f(t,x,p)dp. (34) 

In the case of spherical symmetry Batt [I] gave a global existence proof in 
1977, but later on this was also obtained for general data, independently 
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by Pfaffelmoser and Lions and Perthame ^7] . The method of proof by 
Pfaffelmoser (see also [2H] and ^H] ) is based on a study of the characteristic 
system associated with the Vlasov equation. The quantity that plays a major 
role is Q(t), which was defined in the introduction. If Q(t) can be bounded by 
a continuous function in t, then global existence follows. In jltjj the sharpest 
bound of Q(t) has been given for the system (|33| ) -(|34| ) . where it is proved 
that Q(t) grows at most as C(l+t) log (C + t). As have already been pointed 
out, the continuation criterion in 23 shows that a bound of Q(t) is sufficient 
for global existence also in the case of the spherically symmetric Einstein- 
Vlasov system. It is therefore natural to compare the characteristic systems 
in the two cases, since the natural way to control Q(t) is to obtain bounds 
on the solutions to this system. In the case of the spherically symmetric 
Einstein- Vlasov system this means to bound E(t) = yl + W 2 (t) + F/R 2 (t). 
It follows from the characteristic equations ((18)) and (|19|) that E satisfies 



J o 

A comparison with the analog equation for the Vlasov-Poisson system leads 
to the following observations: 

• the factor (-W) is only present in EV not VP 

• only the m-term is present in VP 

• the dependence on j and p is pointwise 

Having a Gronvall estimate in mind the first point is a severe difficulty for 
the Einstein- Vlasov system since the magnitude of \W\ and E might be of 
the same order, and thus gives much less room for estimating the right-hand 
side. The second and third point concern the regularity of the right hand 
side. The quasi local mass m is an average of the energy density p and 
is thus more regular than p itself, and this fact is important in both the 
Pfaffelmoser- and the Lions-Perthame methods of proof. In the EV case 
there is a pointwise dependence on the matter terms j and p which indeed 
is a difficulty but as will be seen in section 5 this dependence is in a sense 
mild. Other additional difficulties are the following facts: 

• Q 4 instead of Q 3 bound the matter terms p, and p 




(35) 



Here 
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• conservation of ADM mass does not lead to an IP bound of p with 
p > 1 

The densities in the two cases are (in Minkowski space) given by 

p VP = f fdp, and p EV = [ y^T+Wfdp. 
Jrz Jr? 

Again, having a Gronvall estimate in mind it is a disadvantage that p EV 
(and the other matter terms) are bounded in terms of Q instead of Q 3 as 
is the case for p VP '. Conservation of energy in the VP case means that 



\p\ fdpdx, 

is bounded and this in turn implies that p VP (t,-) E L 5 / 3 , whereas con- 
servation of the ADM mass does not qualify p EV into any LP space with 
p > 1. In both methods of proof in the VP case such an a priori bound plays 
a fundamental role. We conclude this section by a comparison with the 
relativistic Vlasov-Poisson system. By starting from the relativistic Vlasov- 
Maxwell system, which describes a charged plasma, and assuming spherical 
symmetry the following system is obtained 

d t f + p-V x f + (3E(t,x)-X7 p f = 0, (36) 
£ = V0, Acj> = p:=f R3 f(t,x, P )d P . (37) 



Here p = p/\/l + p 2 and (3 = 1. This system makes sense also without the 
assumption of spherical symmetry as well as with f3 = — 1 and is called 
the relativistic Vlasov-Poisson system with a repulsive {f3 = 1) or attractive 
((3 = —1) potential. This system is not Lorentz invariant and is not a 
physically correct equation (see [2] for a discussion of this point) but from 
a formalistic point of view the only difference from the VP system is that 
the p in front of the V x /— term has been changed to p, and in this sense 
resembles the EV system. Nevertheless, the mathematical results in the two 
cases are very different. The most relevant results in this context can be 
summarized as follows (we refer to |12| . |14j for details): 

• j3 = 1. A general global existence result is still missing for the RVP 
system but the spherically- and the cylindrically symmetric cases have 
been settled. 

• f3 = — 1. Solutions blow up in finite time for a large class of initial 
data. 
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We see that even in the repulsive case much less is known than for the 
Vlasov-Poisson system, and apparently, in the more relevant attractive case, 
the situation is very different. In conclusion, it is far from clear that the well- 
developed machinery for the classical Vlasov-Poisson system can contribute 
much to the understanding of the spherically symmetric Einstein- Vlasov 
system. Furthermore, the blow-up result for the relativistic Vlasov-Poisson 
system with /3 = — 1 indicates that if the Einstein- Vlasov system prevents 
blow-up in finite time, the method for proving it might be quite subtle. 

5 A mixed L 3+ — and weighted L°° approach 

In this section we make assumptions on the matter quantities and show that 
these are sufficient for bounding Q(t). Since these assumptions are satisfied 
when Q(t) is bounded we obtain an improved continuation criterion. From 
the proof it will be clear that the pointwise dependence on j and p in the 
characteristic equation, discussed in the previous section, is in a sense mild, 
i.e. as long as j and p do not contain too sharp and narrow peaks then their 
contributions can be controlled. 

Assume that there exist constants A and B such that for t £ [0, T[, 

r 2 k(t,r)<A, (38) 
IIp(V)||l3+(b3)<S. (39) 

Here is the ball {x : \x\ < R + + T} and L 3+ means that it is sufficient 
that p £ L ?J+e for any given e > 0. Recall also the definition of k in equation 
l)13[). First note that these assumptions are weaker than the assumption 
that Q(t) is bounded on [0, T[. Indeed, from Q it follows that p < CQ(t) 4 . 
Therefore a bound on Q(t) implies that 

p(t,r)<C, on[0,r[, (40) 

which in view of the fact that the severe difficulties are concentrated at 
r = 0, is significantly stronger than (|38[) (recall k < p). That (|39|) follows 
from (|4fljl is immediate. It is also instructive to see that (|38|) is a natural 
hypothesis in view of the fact that 

171 1 r „ r , 

-<-on[0,T[. (41) 

From the boundary condition A(i, 0) = it also follows that 

m 

* 0, as r — > 0. 

r 
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If we assume that there exists a constant C and a positive number a such 
that 

P(t,r) < — , 



and then estimate m/r we get 



r 



< Cr 2 ' a . (42) 



Thus, if such an estimate for p exists then a < 2, which shows that (|38jl 
is natural. However, j or p might have sharp, and then by (|41|) necessarily 
narrow, peaks so that such an estimate is not valid at all. We have not 
been able to exlude that such peaks exist which is the reason for imposing 
assumption (|H8jl. This assumption still allows for sharp peaks in the matter 
quantity q, which is the part of p which is associated to the tangential 
momenta. The potential difficulty of sharp peaks in j and p is of major 
interest to understand and it might be that this difficulty is a particular 
feature of the Einstein- Vlasov system in Schwarzschild coordinates. We 
have the following theorem. 

Theorem 1 Consider a solution of the Einstein- Vlasov system for regular 
initial data with radial cut-off on its maximal time interval [0, T[ of existence. 
Assume that and i9,9j) hold. Then Q{t) is bounded on [0, T[, and T = oo. 

Proof. The quantity that will be considered is E, and the equation along 
a characteristic (R(t),W(t), F) for E reads 



±E(R(s),W(s),F)= { -^p- 



E. 



(43) 

Note that we often choose to write (— W) since the major case of interest is 
when the characteristic is ingoing, i.e. W < 0. It follows that 



log E(T) =log£(0)+ f ■ 
Jo 



~W)(t) 
E(t) 



dt. 



(44) 

The maximal time interval of existence is [0, T[ and we can assume that 
T < oo. Then there exists a characteristic (R(t), W(t),F) and a Tq e [R_,T] 
such that 

Km R(t) = 0, 

t— >Jo 

since otherwise the finiteness of T would contradict the result in |25| . (It 
is not necessary for our arguments that the limit do exist when Tq = T, 
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but due to the causal character of spacetime it exists.) We consider the 
set of all characteristics with this property and we pick a characteristic 
(R*(t), W*(t), F*) in this set with an associated time Tq such that there is 
no other characteristic in this set with an associated To strictly less than 
T *. Hence we have the following properties of (R*(t),W*(t),F*) and T *, 



lim R*(t) = 0, 



(45) 



and 



f(t, 0, w, F) = for all t < T *. (46) 

From now on we drop the superscript * and consider a characteristic 
(R(t), W(t), F) with the properties above. From the equations Q an d (HH) 
we have 

Hr(t,r) = +47rrp)e 2A , 



and by using @ the time integral in 1)44(1 can be written 



-W(t)) 
E{t) 



H r (t,R(t))-X t (t,R(t)) 



E(t) 



dt. 



(47) 



(48) 



Let us denote the curve (t, R(t)), < t < Tq, by 7. By using that 



dR 



E • 



the integral above can be written as a curve integral 

t^l e C-M+A)(*,,) Mt> r)dr + fc^) e (M-A)( t ,,)^ ( ^ r)dt . (49 ) 
F(t) F(t) 

Let C denote the closed curve 7 + C2 + C3, oriented counter clockwise, where 

C 2 = {(t,0) :T >t>0}, 

and 

C 3 = {(0, r ) : < r < R(0)}. 

Let f2 denote the domain enclosed by C. We get by Green's formula in the 
plane 



c 



J5(t) 

~W(t)) 
E(t) 



fj, r dt 



Or 



At 



dtdr 



(50) 
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Note that the factor W(t)/E(t), which really only is defined along the char- 
acteristic, has been considered as a function of t alone, i.e. for given t 
it is defined to be constant in r with the value it has at (t,R(t)). Alter- 
natively, E could have been viewed as a function of t and r, by taking 
E = -\Jl + W 2 (t) + F/r 2 . The final conclusion would of course be the same. 
An interesting fact is now that the first term in the right hand side which 
contains a combination of second order derivatives of [x and A can be sub- 
stituted by one of the Einstein equations ©• Indeed, we get 



d t ( e-^ +A A f 



dtdr 



-W(t)) 
E(t) 



(-wit)) 

E{t) 
-W(t)) 



X r ))e^ x - X tt - (m - At)) e-^ +A 



dtdr 



47r(/e 



(fi r — A r )- 



,/^+A 



E(t) 



4ir(p — p) + Airq 



2m 



dtdr 
dtdr. 



(51) 



Here we used (|47|) for ^ r and that 

A r = (4-irp 



~k)g , 



which follows from ©. By using the characteristic equations (j!8|) and ()19|) 
we get that the second term in (p)U|) is given by 



E 2 



Fe^~ 



R 3 E 2 



Xtdtdr. 



(52) 



To summarize we have obtained the following identity 



C 



~W(t)) 
E(t) 



e-» +x \ t dr + 



E(t) 1 



,M~A 



/i r dt 



//J(— 



=:T! + T 2 . 



(53) 



Note that the function in square brackets in T 2 only depends on t and its 
argument is (t,R(t)), whereas it is (t, r) for e~^ +A At. The term T\ must 
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be considered as the main term. If the integrand in (|48[) is nonnegative 
we can drop the factor (— W(t))/E(t) and estimate the time integral by T\ 
alone (as in section 7). More importantly, T\ has a geometrical meaning; 
ignoring the factor {—W)/E it is the integral of the Gauss curvature of the 
two dimensional manifold Rf x with metric ds 2 = —e 2fl dt 2 + e 2X dr 2 . (We 
really have in mind the formula in (|51|) obtained before we used Einstein's 
equations to get an expression in terms of matter quantities.) It is interesting 
to note that the principal term in the work on BV solutions of the ESF 
system by Christodoulou 7\ is of the same form, and in complete analogy 
with our result (see p. 1172 in [7j). 

By our assumption that p £ L 3+ , it is an easy matter to estimate T\ as 
it stands. However, we may first note that in the case of initial data with 
angular momentum cut-off, T\ can be simplified since the following relation 
holds 

/■oo /-oo n 

p(t,r)-p(t,r)=q(t,r)+ / -j^f(t,r,w,F)dFdw, 

where the last term is harmless. Indeed, note that 

/oo poo 
/ r, w, F)dFdw < —r 2 p, 

~ oo J 

where we used the property that F > F_ . This term is bounded in terms of 
the ADM mass. Therefore, the remaining and the principal part, Tf Tl , of 

' j L p g Qjd. s 



<> 



E(t) V r3 J 



which has the advantage that only q and not p is explicitly in the integrand 
(m depends of course on p). This will have significance in section 7. 

We now estimate T\ in ()53|) . By using \ W\/E < 1 and e^ +x < 1, we have 

|2i| < J J + 8irp\ dtdr = + 8irpJ drdt. 

(54) 

It should be noted that it is our assumption ()39|) on p which allows for this 
rough estimate. In a more careful treatment the major difficulty should be 
associated with an ingoing trajectoria, i.e. W < 0, and therefore the signs 
of the terms in the integrand should be important (see section 7). By our 
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assumption (|3*§)) we get by Holder's inequality (1/p + 1/p' = 1, p > 3) 

/•r / /-r \ i/p' / rr \ i/p 

m = rf p{t,rf) dn < I / r/ 2 <i?7 I I / rf^dn) 



i/p 

< Cr 2+ ( / rffPdrj ) , (55) 



and 



R(t) ( f R(t) \ 1/P ' f fR(t) \ 1/P 

p(t, r)dr<\ j r~ 2p ' p dr \ f J r 2 ffdr \ 

R(t) \ l 'P 



o 



<CyJ rYdrj , (56) 

where we have used that p'/p < 1/2. In view of these estimates we obtain 

|Ti| < CBT. (57) 

Now we consider the term T2. We begin with a lemma showing that our 
hypotheses imply that A is bounded. 

Lemma 4 Assume that J3SJ) and \'jy\) hold on [0, T[. Then A is uniformly 
bounded on [0, T[. 

Proof. From the relation (1221) we see that A is bounded as long as m/r is 
strictly less than 1/2. From (|39|) it follows that 

m r 11 1 r ^ 

— < ^tjIIpIIl 3 < ^7J- 

If r < 3 2 / 3 /4i? it follows that m/r < 1/4 and A is bounded by log 2/2. For 
r > 3 2 / 3 /4i? we consider equation (jHJ) and in view of (|3*5)) we get 

Air A 16vrAB 

and since A is bounded initially it follows that A is bounded on [0, T[ when 
r > 3 2 / 3 /45. By combining these estimates the lemma follows. 

□ 
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To estimate T2 , note again that the argument of the part of the integrand in 
square brackets is (t, R(t)) while the argument of e~ M+A At is (i, r). By using 

H r (t,r) = (m/r 2 + 47rrp)e 2X , 

we obtain four terms in the integrand of T^. Since 

\ t (t,r) = -4vrre M+A j, 

and the sign of j is indefinite we cannot drop any of these terms. For the 
first term we have 

e (-M+A)(t,r) e (M+A)(t,r) 47rT (_ J .( tj r )) dtdr 

e 2X Air\j\ dtdr. (58) 

Here we used that r < R(t). The expression in square brackets is bounded 
by 1/2 and since A is bounded in view of Lemma 5, this term is bounded by 

C I I \j\drdt < C||p|| X 3+T < CBT. 
J Jn 

The remaining terms are estimated in a similar fashion by again using that 
r < R(t) together with (|38]> to estimate the terms rR(t)p(t, R(t)), and 
rR(t)\j\. In conclusion we have obtained that 

Ti + T 2 < CBT. 

Therefore, a bound on E{t) follows from (|44jl if we can bound the curve 
integrals over C2 and C3. The curve integral over C3 is obviously bounded 
in terms of the initial data and the curve integral over C2 is zero by (|46j) . 
This completes the proof of Theorem 1. 

□ 

Remark. We have seen that our assumption on p makes the estimate easy 
and we have only to consider the r— integration in the T\ term. It is however 
clear that a more careful treatment must benefit from the time integration. 
Indeed, if p behaves like r~ 3 / 2 , then an approach that only focus on the 
r— integration in the T\ term will fail while such a behaviour of p is perfectly 
allowed by (|23|) and the a priori bound in Lemma 1. The methods of proof 
of global existence in the case of the Vlasov-Poisson system make use of such 
a time integration in a cruscial way. 



m(t,R(t)) 

. R(t) 2 
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6 An approach along the characteristics 



In this section we give conditions on /j, r and A r along a characteristic which 
guarantee the boundedness of Q(t). These conditions constitute a borderline; 
we will show that when these conditions are satisfied, W(t) must become 
before R(t) = 0, and thus the characteristic "turns" from ingoing to out- 
going, and if they are not satisfied, we show under additional hypotheses 
that e^ _A becomes sufficiently small to make R(t) stay uniformly away from 
zero. These results clearly support the idea that there is a mechanism which 
prevents blow-up of the solutions. 
Define 

ff ( r ) = I(l + _L) ) <r < 1/3. 
r logr 

We make the following hypotheses. Along a characteristic (t, R(t), W(t),F) 
we assume that the following conditions hold when R(t) < 1/3, 

lir(t,R(t))<g(R(t)) (59) 

fr(t,R(t))-)*(t,R(t))<j^. (60) 

The condition R < 1/3 is only technical and we could have taken any R < C 
by minor changes of the presentation. We are going to show that Q(t) is 
bounded under these hypotheses but first let us discuss these assumptions 
and put them into a context. In spirit these assumptions are similar to the 
hypothesis Q38JI in the previous section. Observe that 

= (JS + 47rrp)e 2A , 
and we see that 1)59(1 implies that 

Thus, since g(r) is roughly 1/r this is of the same order as 1)38(1 with the 
main differencies that the constant A is now precisely specified and that the 
factor e 2A is included. In this situation the factor e 2A is actually bounded 
by 3, since by (JSHJ) (for R < 1/3) 

and using fl2D it follows that m/R < 1/3 so that e 2A < 3. 

A natural question is what can happen if these assumptions are not 
satisfied and the following lemma gives a partial answer. 
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Lemma 5 Let (R(t),W(t),F) be any solution to the characteristic system 
on a time interval I = [ti,tz] with 1/3 > R(t{) > 5 > and W(t) < on I, 
and assume either that 

fi r (t,R(t)) > g(R(t)) and fi t (t,R(t)) < 0, 

or that 

fi r (t, R(t)) - X r (t, R(t)) > -— r and fi t (t, R(t)) - X t (t, R(t)) < 0, 

H(t) 

on I. Then there exists an e > 0, which only depends on 5 such that 

R(s) > e, s e [h,t 2 ]. 

It should be noted that by making assumptions on the signs of both fit and 
At along the entire characteristic (note that in Lemma 5 these assumptions 
are only made on that part where the spatial derivatives are large) one can 
easily get control of ingoing characteristics with F > 0. We formulate this 
in a lemma. 

Lemma 6 Let (R(t),W(t),F) be any solution to the characteristic system 
on a time interval I = [ii,t2] with R(t\) > 5 > 0, F > 0, and W(t) < on 
I. Assume that fj,t(t,R(t)) < and that Xt(t,R(t)) > 0. Then there exists 
an e > 0, which only depends on 5 such that 

R{s) > e, s e [h,t 2 ]. 

Before giving the proofs of these lemmas let us discuss the sign assumptions 
in Lemma 5 and compare the content of these two lemmas. The assumption 
that fit < when fi r > g (or fit — X t < when fi r — X r > 1/R) is of course 
the weak part of lemma 5. However, note that it always holds that fi < (or 
fi — X < 0) and that for collapsing matter where the density is increasing, 
fi (or fi — A) would decrease and in this sense the case fit{t,R(t)) < 
(or fjLt(t,R(t)) — X(t,R(t)) < 0) should correspond to a regime where most 
matter collapse and little disperse. At least intuitively, this should be the 
worst case. Since the sign assumption in Lemma 5 is only made on the part 
of the characteristic where the spatial derivative is large it strengthens this 
picture since it is in the collapsing regime that these derivatives are expected 
to be large. This is indeed supported by the numerical simulation in |26j . 
Moreover, in Lemma 6 we have taken sign assumptions on both fit and At 
instead of only on fit (or on fit — At). Lemma 6 also requires initial data with 
cut-off (F > 0), whereas the proof of Lemma 5 holds generally. From this 
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discussion we therefore consider Lemma 5 as being of more interest than 
Lemma 6. 

Proof of Lemma 5. Let us first consider the case when [i T > g(R). 
Writing the characteristic equation for R(s) as 



dRr\s) 
ds 



(-W(s)) 



R(s) E(R,W,F)_ 



R-\s), 



and noting that 

{-W)/E < 1, 

we see that can be controlled if we can show that 

(v-WMt)) i 
< Clog 



R(t) 



R(tY 



Now 



^- x ){s,R{s)) e^(s,R(s)) 



R(s) 



< 



R(s) 



since A > 0, and we have 
ds~ Rjs) 



dR\ e» fWe^~ x \ 

~R 



< 



Hr(t,R(t)) 



1 



1 



R(t) 



ER 

ER 

dR.e^ 



R(t) log (l/R(t)) K ds' R' 
Here the assumptions on \i r and \it were used. Hence we get for t € [ti,t2], 



e n(t,R(t)) e /i(ti,fi(ti)) l ogj R(t) 



R(t) 



R{h) logfl(ti)' 



which completes the first part of the lemma. In the second case the proof is 
analogous but easier since it follows that 



so that 



d / e n(s,R(s))-Ms,R(s)) 
~ds~ y Rjs) 

e (n-\)(t,R(t)) 



< o, 



R(t) 



is nonincreasing on [ti , *2] • 
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□ 



Proof of Lemma 6. The spirit is the same as in the previous proof. Since 
F > we have 



\W(s)\ 



< 



< 



R{s) ^/\ + W 2 + F/R 2 



R-Hs) 



F 



E(R,W,F)R- 1 (s) 



This inequality shows that it is sufficient to bound e^E. By using the char- 
acteristic equation for E in (|43j) and inserting the expressions for fi r and At 
and using the assumptions of the lemma it is an immediate fact that 

d 



ds 



(e 



E) < 0. 



□ 



We now state and prove the main result of this section. 

Theorem 2 Let f E Ir_,F- and let (R(s, 0, r, w, F),W(s, 0, r, w, F),F) 
be any solution to the characteristic system U8\) - U9\) on [0,T[ such that 
fo(r,w,F) ^ 0, and assume that \59\) and \6fl\) hold along this characteris- 
tic. Then Q{t) is bounded on [0, T[. 

Proof. We can without loss of generality assume that F_ > 1 . The argument 
gets more transparent with this assumption but it is clear from below that 
it is easy to modify the argument for any F_ > 0. Let to = T — 1/6. Since 
the solution is regular on [0, T[ we have that Co := Q(io) < °o> in particular 

R(t ) > and \W\ < C , for all characteristics (R(t),W(t), F). We 

now consider the time interval [io)^"[ an d we show that all matter stays 
uniformly away from the centre on [io)^[- The class of characteristics with 
R(t ) > 1/4 will have R(t) > 1/12 for all t G [t ,T[ by the fact that 
\dR/ds\ < 1. Hence we consider the complementary class of characteristics 

with ^— < R(t ) < 1/4. First we note that if dR/ds > then dR/ds / 0, 
i.e. W ^ 0, as long as R < 1/3 since we have for any t with W{t) = and 
R(t) < 1/3, 



e fl - x n r (l + F/R 2 ) + 



> 



dW 
ds 

-1 

~R + R\og(l/R) 



Fe^ 



-x 



1 



+ 



F 



R 3 log(l/i?) 



R 3 E 

> 0, 
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by (|59|) and since F > 1. Thus by continuity, an outgoing characteristic 
cannot turn around and become ingoing as long as R < 1/3. If it turns back 
when R > 1/3 it can never get closer than R = 1/4 for t <T. We have thus 



d 
ds 



found that the case of interest is when dR/ds < and < R(to) < 1/4. 
We now show that for any such characteristic, W will become zero before R 
becomes zero. By (|T9*)) we get 

±{e 2X W 2 ) = -2\ r —e»+ x W 2 - + W 2 + F/R 2 ) + 2WF f +X 

ds EE R 6 E 

x sdR o\ TT .o IF , 1 .dR n\ n dR 9 \ 
= -2(/x r - \ r )—e 2X W 2 - -^(pr - ^)^e 2A - 2^e 2 Vr- 

Hence, by (|59|) . and the fact that logR < — 1 and F > 1 we have 

R^ r R } ds 6 ds 6 llr ) e 

< ( 7 2F - 2» r ) e^—efo^-^ds 
~ \R 3 logR P J ds 

~ R 3 log R ds 

-F dR 2 -F dR 
~ R 3 log R~ds~ ~ RlogRls~' ^ 

Integrating this inequality gives 

< -F log | log fl(t) | + F log | log i2(t ) I , (62) 

so that 

e 2X(t,R(t)) W 2^ < 
< e -/4^-Ar)#* ( e 2A( f0 ,^ ))^ 2( } _ F1 | logfi(t)| \ 

This relation implies that W will vanish before becomes 0. From the 
arguments of the first part of the proof we have seen that dW/ds > when 
W = 0, so that the characteristic necessarily becomes outgoing and cannot 
turn back in again as long as R < 1/3. 

□ 
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7 The case of purely ingoing matter 



A special case of interest is when all matter is falling inwards. This is for 
example the situation when dust is used as matter model; if the matter 
is initially falling inwards, there will be no outgoing matter at any time. 
Moreover, in the numerical simulation j^f)] most matter is seen to be falling 
inwards in the supercritical case when black holes (presumably) form. We 
believe that a good understanding of this case will also have impact on the 
general case since control of ingoing matter give (in a certain sense) control 
of the outgoing part. In this section we will take initial data such that 
all particles are initially falling inwards and we will assume that all matter 
continue to fall inwards at least on some time interval [0, T[. The other 
assumption is connected to the relation between the enclosed mass m(t,r) 
and the areal radius r. It is well-known that there are closed null geodesies 
when r = 3M in the Schwarzschild spacetime with mass M. Therefore it 
is not surprising that the relation r = 3m turns out to have significance 
also in the Einstein- Vlasov system. We will show that when all matter is 
falling inwards, control of ingoing characteristics can be obtained either if 
r < 3m(t, r) everywhere, or if r > 3m along the characteristics. Before 
stating our hypotheses we define the outgoing and ingoing part of j, 

j + (t,r) 

On the time interval [0, T[ we assume that 

j+ = 0, (66) 
r > 3m(t, r). (67) 

We remark that the first hypothesis can be replaced by assuming that j + 
is bounded but we find it more transparent to think of the results in this 
section as the case of purely ingoing matter. The first question we ask is 
what can happen if the second hypothesis is not satisfied and the following 
lemma gives a partial answer. 

Lemma 7 Let (R(s),W(s),F) be any solution to the characteristic system 
on a time interval I = [ti^] such that R(t\) > 5 > 0. Assume that 166]) 
holds and that 

R(s) < 3m(s,R{s)), 



7T 
7T 



dw / dF wf(t,r,w, F), 
Jo 

/0 poo 
dw / dF wf(t,r,w, F). 
-oo JO 



(64) 
(65) 
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on I. Then there exists an e > 0, which only depends on 5 such that 



R(s) > e, s e [h,t 2 \- 
Proof. The proof is similar to the proof of Lemma 5. Indeed, recall that 

f°° m(t n) 
r) = - e 2A ( v \ " + 4nrrip(t, r,))dr]. (68) 

Jr V 

Let us write 

H = p, + p,, (69) 

where 

H{t,r) = -J e —^—di). 

Since e' 1 < we proceed as in Lemma 5 and consider -^(e^/R). From the 
equations © and (j221) it follows that 

Mt,r) = r4nje x+ ^l + 2e 2X ^^)dr ] , (70) 

Jr V 

so /it < by (|66|) . It is then easy to show that e^/R is nonincreasing by a 
similar calculation as in Lemma 5. 

□ 

The main result in this section concerns an ingoing characteristic for which 
d(e^E)/ds > 0. We have seen above that control of this quantity is sufficient 
for controlling R(t) . The condition that the time derivative is nonnegative 
should then be the relevant case. One can think of dividing the time interval 
into subintervals where e^E does increase and where it decreases and then 
add up the contributions where d(e^E)/ds > 0. However, if the sign of 
d(e^E)/ds changes an infinite number of times on [0, T[, then it is not clear 
that the infinite sum of the contributions converge due to the nature of the 
estimate we obtain below. This is the reason for including this assumption. 



Theorem 3 Let /q 6 and such that fo(-,w,-) = if w > 0. Let 

(R(s, 0, r, w, F), W(s, 0, r, w, F),F) be any solution to the characteristic sys- 
tem J72P-J23) on [0,T[ such that f (r,w,F) ^ 0. Assume that d{e (l E)/ds > 
on [0, T[, and assume that \66\) and \61\j hold. Then there exists an e > 
such that 

R{s, 0, r, w, F) > e, for all s £ [0, T{. (71) 
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Proof. We will only sketch the proof since it is to some extent similar to the 
proof of Theorem 1. Again we consider e^E and we show that this quantity 
can be controlled. We have 

— (E(R(s),W(s),F)e^ s ' R ^ 
as V 



47T 



+A 



W 2 



\w\ 



R(s) 



ef+lj ( i + ^1 ) drj 



Ee^. 



Here we have inserted the formula (|7U|) for Note also that the terms 
involving m/R? have cancelled. From the assumption that e^E is nonde- 
creasing we have (note that the integral term is nonpositive by (|66j0 

W 2 \W\ 



and thus 



We get 



E 2J+ E 



W\ \W\ 

P < —^3 + P- 



E 



^(E{R{s),W{s),F)e^ R{ - s)) 
as \ 



< 4vr 



Re? 



+A 



\W\ 
E 



J+P + 



e^ +X j{s,i!) 1 + ) d V 



R(s) 



>1 



Ee». 
(72) 



This is a similar to equation (|43j) and we are again going to apply Green's 
formula in the (t,r)— plane. This time we will choose the outer domain f2, 
which is enclosed by the curve (t,R(t)), t G [0, T[, which we denote by 7, 
together with the curves 

C 2 = {(T,r):R(T)<r<R°°}, 
C 3 = {(t,R°°) :T>t>0}, 
C 4 = {(0,r) : R°° > r > R(0)}. 

Here R°° > R+ + T, so that / = when r> R°°. The closed curve 7 + C 2 + 
C3 + C4, we denote by C and it is oriented clockwise. We integrate (f?2|) in 
time and obtain a curve integral over 7, analogous to ()49|) . and we note that 
the part containing \ W\j/E +p can be written 



X t dr + e M x jl r ds, 
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and we then have by Green's formula 

j> e-^ +x \ t dr + e^ x jl r ds 
d t (e-» +x \ t ) - d r (e^fir) dsdr 



d t (e-^Xt) - d r (e^~V) dsdr + I I d r (e^ +A ^) dsdr 



(73) 



Using the calculations that led to equation (|51[) we obtain the identity 

-^ +x \ t dr + e^ x fi r ds 



'C 



I'll 



(p + p)e 2X — + 2p- p 
r 



drds 



+ 



ft J-oo JO 



r 2 E 



f(t,r,vu, F)dFdw drds. 



(74) 



Recall from (|72jl that we wish to compute the curve integral 

r T roo 



(I v E 



\ pJ. poo 

+ p) + / ^je x+t "{l + 2e 2X —)dr]ds 
/ JO JR(s) V 



lR(s) 
rT roo 



I e-^ +x \ t dr + e»- x fi T ds + [ [°° 47tfe A+ "(l + 2e 2X -)dr,ds 
Jc Jo Jr(s) V 



la. 



Xtdr + e' 1 x fi r ds. 



'C2+C3+C4 
Since j = when r > R°° we get 



(75) 



/* f'T f'OO 

<b e-^ x X t dr + e^- x jl r ds+ / / 4irje x+tl (l + 2e 2A 



47re M 



+A 



( p + p)e 2X ^ + 2p-p + j(l + 2e 2X ^) 



-)drjds 



drds 



n J-00 Jo 



00 47r 2 e/ ,+A 
r 2 £ 



/(i, r, u;, F)dFdw drds. 



(76) 



The definitions of the matter terms and assumption Q66J) imply that 

P < "i < P, (77) 
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and the first term in (|76|) can now be estimated by 

(p + p)e 2X — + 2p- p + j{\ + 2e 2X — ) 
r r 

= ~{P ~ P)(l " e 2A -) + (p + + 2e 2A -) < 0, 
r r 

where we used 

9 N,m m/r 

e — = - — r < 1) 

r 1 — 2m/r 

by assumption (|67|). A bound for the second term in (|76jl follows as in the 
proof of Theorem 1 where Tf nn was introduced. The two first terms in (|75j) 
have now been estimated and we are left with the boundary term 

- / e'^Xtdr + e fl ~ x il r ds. 

JC2+C3+C4 

The integral over C2 is nonpositive by assumption since 

r r-R 00 rR°° 

- / e-^ +x X t dr + e fl - x p, r ds = - e"^ +A A t dr = Aite 2X rj < 0. 
Jc 2 JR(T) Jr{t) 

The integral over C3 is zero by letting R°° tend to infinity since m < M. 
Finally, the integral over C4 depends only on the initial data and is thus 
bounded (we remark that it is at this point it is troublesome to add up an 
infinite sum of contributions as was previuosly discussed) which completes 
the proof of the theorem. 

□ 



Remark. An alternative way of presenting a similar result would have been 
to assume that the sign of d(e >l E)/ds, and of R(t)—3m(t, R(t)) only change a 
finite number of times on [0, T[. These assumptions together with (|66JI would 
be sufficient for global existence. Such a proof would rest on the results in 
this section together with a more involved analysis of the boundary terms. 
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